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ABSTRACT 

Singlet Higgs bosons present in extensions of the MSSM can have sizable Yukawa couplings 
to the b quark and the r lepton for large values of tan (3 at the 1-loop level. We present an 
effective Lagrangian which incorporates these tan /3-enhanced Yukawa couplings and which 
enables us to study their effect on singlet Higgs-boson phenomenology within the context 
of both the mnSSM and the NMSSM. In particular, we find that the loop-induced coupling 
can be a significant effect for the singlet pseudoscalar, and may dominate its decay modes. 
Further implications of the tan /3-enhanced Yukawa couplings for the phenomenology of the 
singlet Higgs bosons are briefly discussed. 
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1 Introduction 



The Minimal Supersymmetric extension of the Standard Model (MSSM) provides a self- 
consistent framework to technically address the gauge-hierarchy problem related to the 
Standard Model (SM) Higgs boson. The 2-Higgs doublet potential of the MSSM is highly 
constrained at the tree level which makes it very predictive. Radiative corrections from 
third generation squarks affect significantly the mass of the lightest neutral Higgs boson 
Hi and generically shift it upwards by an amount of ~ 30-40 GeV. The Hi boson has to 
be lighter than about 135 GeV for large values of tan/3 > 20, where tan/3 is the ratio of 
the two Higgs vacuum expectation values (VEVs). For low values of tan (3 ~ 3, most of the 
parameter space that gives a mass for the lightest neutral Higgs boson above the present 
LEP experimental limits has been excluded. For a comprehensive analysis performed by 
the LEP Higgs working group, see [1]. 

One theoretical weakness of the MSSM is the so-called /i-problem [2,3]. This is related 
to the fact that the /i-parameter describing the mixing of the 2 Higgs superfields in the 
superpotential, i.e. {iH u H d , has to be of order the soft SUSY-breaking scale M SU sy ~ 1 TeV, 
for successful electroweak symmetry breaking. Instead, within the context of supergravity 
(SUGRA), the /i-parameter is in general not protected by gravity effects and so expected 
to be of order Planck mass Mp\. A natural solution to the /x-problem may be obtained in 
extensions of the MSSM, where the /i-term has been promoted to a dynamical variable, 
e.g. to a SM-singlet chiral superfield S. In such a setting, the scalar component S of S 
generically acquires a VEV of order Msusy, thereby giving rise to a /x-term of the required 
order. 

Since the /i-term is replaced by the term \SH u Hd in singlet extensions of the MSSM, 
the resulting superpotential will exhibit an unwanted global Peccei-Quinn (PQ) symme- 
try [/(1)pq, unless further additions or assumptions are made to the model. The PQ 
symmetry must be explicitly broken to avoid the appearance of visible electroweak-scale 
axions after the spontaneous symmetry breaking SSB of £/(1)pq. The choice of discrete 
and gauged symmetries used to break the PQ symmetry distinguish between the different 
models which have been studied in the literature [3], including the Next-to-Minimal Super- 
symmetric SM (NMSSM) [4], the minimal nonminimal Supersymmetric SM (mnSSM) [5], 
the £7(l)'-extended Supersymmetric SM (UMSSM) [6] and the secluded L r (l)'-extended 
Supersymmetric SM (sMSSM) [7]. 

Irrespectively of the details of the particular model, singlet Higgs bosons have no tree 
level couplings to SM fermions or gauge bosons. It has long been known [8, 9] that within 
the MSSM, threshold corrections to the Yukawa couplings to b quarks and r leptons can 
become significant in the limit of large tan f3, which partially overcome the loop suppression 
factor. In regions where the mixing between the Higgs particles is negligible, the one-loop 
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correction can dominate the Hi — > bb decay width [10]. 

In this paper we show that an analogous tan f3 enhancement takes place for the Yukawa 
couplings of the singlet Higgs bosons in minimal extensions of the MSSM. In particular, we 
explicitly demonstrate how the effective couplings Sbb and St + t~ are generated radiatively 
through squark-gaugino loops and their size can be significant, e.g. of order the SM Yukawa 
couplings. In the limit in which the doublet decouples from the low energy spectrum, 
these one-loop couplings provide the dominant decay mechanism for light singlets. Recent 
work has proposed the possibility of light Higgs singlets as a solution to the little hierarchy 
problem [11]. In this scenario the threshold corrections are not only important for the 
Higgs searches themselves, but they may potentially provide a first signpost towards the 
physically realized region of SUSY parameter space. 

The paper is organized as follows: in section 2 we present the effective Lagrangian for 
the Higgs-boson couplings to bb and to r + r~, along with analytic expressions for the dom- 
inant contributions. Section 3 discusses the phenomenological implications of the radiative 
singlet-Higgs Yukawa couplings for the mnSSM and NMSSM. Our notations and conven- 
tions regarding the Higgs sector follow those of the first paper in [5], whilst those regarding 
the squarks and charginos are outlined in Appendix A. Our conclusions and possible future 
directions are presented in Section 4. 

2 Effective Yukawa Couplings 

In this section, we derive the general form of the effective Lagrangian for the self-energy 
transition fi — > fa, in CP-conserving Higgs singlet extensions of the MSSM, where / may 
represent a 6-quark or a r-lepton. We then use the Higgs-boson Low Energy Theorem 
(HLET) [12, 13] to compute the 1-loop effective Higgs-boson couplings to b quarks and r 
leptons. 

The general effective Lagrangian for the self-energy transition fi — > fn in the back- 
ground of non-vanishing Higgs fields may be written down as 

- CL = h f f R ^T + A/pS, + h.c, (2.1) 

where $° 2 = 75 ( V J.,2 + $1,2 + ^1,2) are the electrically neutral components of the two 
Higgs doublets -ff^Jj, S = (vs + 4>s + i> a s) is the singlet Higgs field and the functional 
A/[$5j S] encodes the radiative corrections. Given that the VEV of the effective La- 
grangian — £{ elf should equal the fermion mass rrif, an expression for the effective Yukawa 

1 Here we adopt the convention for the Higgs doublets: H u = <£> 2 , Hd = iT2&*, where r 2 is the usual 
Pauli matrix. 
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coupling hf can be found, i.e. 
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where (Af) is the VEV of Af, where the renormalization scale is set at Mgugy. The 
contributions to the functional Af which get enhanced at large tp = tan/3 have been 
well understood within the framework of the MSSM [9, 10]. As is shown in Fig (TJ the 
dominant contributions to Af come from diagrams with gluinos and bottom squarks and 
with chargino and top squarks in the loop. 

We may relate the self-energy effective Lagrangian C[ cl{ to the effective Lagrangian for 
the Higgs-boson couplings to the fermion /, using the HLET [13]. In terms of the physical 
Higgs fields -£^1,2,3 and A\^, the effective interaction Lagrangian reads: 
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Here the orthogonal matrix H (0 A ) is related to the mixing of the CP-even (CP-odd) 
scalars and the loop corrections are given by the HLET 
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In (12 .4p and (12. 5p . we have neglected the one loop contributions to the <pi coupling, since 
they are small, i.e. A^ 1 < Af' s tp. 



2.1 Effective 6-quark Yukawa Couplings 



As in the MSSM, there are tg-enhanced contributions to the self-energy of the b quark 
from both gluino and chargino exchange diagrams, as shown in Fig. [TJ Evaluating these 
^-enhanced diagrams at zero external momentum and neglecting subdominant terms pro- 
portional to a w yields 
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Figure 1: Dominant contributions to the bottom quark self energy at large tp, in the single 
VEV insertion approximation. 



In the above, I(a, b, c) is the usual 1-loop integral function given by 

ab\n(a/b) + bchxib/c) + acinic/ a) 

I{a,b,c) = — r , 2.8 

(a — b){b — c)[a — c) 

and V, U are the chargino- mixing matrices defined in Appendix A. Note that V, U are 
functionals of $5 2 and S, as are the sbottom quark masses , stop quark masses 2 
and chargino masses 2 . Explicit expressions for the masses and mixing angles are given 
in Appendix A. 

We may then use the HLET to calculate the corresponding graphs with an addi- 
tional zero momentum Higgs insertion. The lowest order <psbb graphs are shown in Fig. [2j 
Neglecting again terms proportional to a WJ the coupling parameters Af 2 ' s are given by 
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where A b = — ^utg and A 4 = A 4 — njtp. The derivatives act on all Higgs-dependent 
functionals to their right, generating a rather lengthy expression which we do not show 
here explicitly. These derivative terms represent higher number of Higgs insertions, beyond 
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Figure 2: Dominant contributions to the couplings <f)sbb, 4>s r+r , a sbb and a s r + T , in the 
large tp limit. 



the usual single Higgs insertion approximation often followed in the literature. They can 
be t/3-enhanced in certain regions of the parameter space, especially for the chargino case, 
and are therefore consistently included in our numerical analysis in Section 3. 

The presence of the singlet in the model does not alter the 1-loop tan j3 enhanced 
couplings of the doublet Higgs fields well-known from the MSSM [14]. As a consistency 
check, we have compared the gluino-exchange terms of (12.91) with corresponding results from 
the first reference of [14], and we find that the HLET calculation gives identical results. 

Correspondingly, the loop-induced coupling parameters A^ 2,s are given by 
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Notice that both $° an d the singlet field S enter the 6-squark masses only through the 
mixing term (Ab<&± — A5$o)- As a result of this, it is easy to derive that dm\ /dips — 
(1)2/ vs) dmf /d(f>2- With the aid of the latter, we then find that At 2 is related to At s by 
A^ s = ^| Af 2 , if the gluino-exchange diagram is only considered. In the same approxi- 
mation, a similar relation holds true for the pseudoscalar 1-loop couplings: A^ s = ^ A% 2 . 
However, this simple scaling behaviour is broken by the chargino-exchange diagram, as 
the background Higgs fields $° an d S enter the chargino masses and mixing angles in a 
non-linear manner [cf. (1A.12I) and (1A.13I) . 

Finally, it is worth commenting on the fact that the coupling of the neutral would-be 
Goldstone boson G° to the b quark is proportional to the 6-quark mass, as a consequence 
of a Ward identity involving the Zbb-coup\ing. Specifically, the coupling parameter Af°, 
which is computed by 

Af = ,V2 (^A t ) 

is given by 

Consequently, the G°66-coupling has the tree-level SM form in the limit of zero external 
momentum. 



(2.14) 



2.2 Effective r-lepton Yukawa Couplings 

The derivation of effective r-lepton Yukawa couplings goes along the lines discussed above 
for the 6-quark case. At the one loop order, there is now only one t/j-enhanced diagram con- 
tributing to A T , which originates from the chargino-Higgsino-exchange diagram of Fig. [2j 
The effective functional A T pertinent to the r-lepton self-energy is given by 
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where Mjr is the soft SUSY-breaking mass term for the left-handed sleptons and the loop 
function B°(p 2 , a, b) at p 2 = takes on the simple form 
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where Q is the renormalization scale, which is conveniently taken to be Q' — ul^ 2 ±v± l . 

Notice that the non-holomorphic couplings both receive renormalization scale-dependent 
contributions. This should not be suprising, as the couplings also contain contributions from 
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the mixing of 02, s with the holomorphic <j}\. We have checked that all scale-dependent terms 
vanish in the zero mixing limit of <p2,s —> 4>i, corresponding to V\ — > 0. 

As was done above for the fr-quark case, the 1-loop Higgs-boson couplings to t + t~ 
may be computed from (12.151) by means of the HLET, where the dominant 0st + t~ dia- 
gram is shown in Fig. [2j In extensions that include right handed (s)neutrinos, a second 
diagram mediated by Higgsino exchange must also be considered. This new contribution 
may be significant if (s)neutrinos are not too heavy. Such further extensions may be studied 
elsewhere. 



3 Phenomenological Discussion 

In this section we analyze the implications of the loop-induced parameters A^ s and Aj s 
for the Higgs-boson couplings and for the Higgs-boson phenomenology in general. As was 
already mentioned, the 1-loop coupling of the singlet Higgs boson to the b quark and the r 
lepton becomes significant at large values of tan/5. In our analysis, we adopt a benchmark 
scenario where the singlet Higgs-boson effects get enhanced. Unless is stated otherwise, 
the default values of the SUSY parameters for our benchmark scenario are 

tp = 50, 
M z = 90 GeV, 

M~ b = 110 GeV, M f = 200 GeV, (3.1) 

A b = 1 TeV, A T = 1 TeV, 

M 2 = 600 GeV, M 3 = 400 GeV . 

Notice that an important constraint on the choice of the above parameters comes from the 
LEP data. This constraint is included in our analysis. 

Given the model parameters (13.11) . the coupling-parameter ratios Ay S /Aj 2 and Aj s /A' 
are shown in Fig.[3j as functions of the supersymmetric coupling A, keeping the /x-parameter 
fixed. Since the radiative corrections to the Yukawa couplings are dominated by SUSY 
QCD effects in the case of the quarks, one might expect for the ratios to be approxi- 
mately given by Af s /Af 2 pa A^ s /A^ 2 pa v/vs- Specifically, this ratio should reach the 
value 1 for A ~ 0.65. As is illustrated in Fig. [3l whilst this approximation is valid for the 
psuedoscalar couplings, the subdominant corrections to the scalar couplings do not share 
this simple scaling behaviour and so give rise to a somewhat different relative magnitude 
for Af s /Af 2 . Clearly, including interference effects between the contributing terms, the 
coupling parameter A^ s becomes comparable with Af 2 only for large values of A. 



H = -j^Xvs = 110 GeV, 

= 300 GeV, 
Mf = 600 GeV, 
A t = ITeV, 
Mi = 400 GeV, 
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Figure 3: Numerical estimates of the ratios A^/A^ and A^/A^ 2 T as functions of X. The 
values of the model parameters are given in (COP- 



As can be seen from the effective Lagrangian f)2.3p . the physical Higgs-boson couplings 
to the b quark and the r lepton, e.g. Hi 2 ^ff and Ai ^ff, consist of two contributions. The 
first contribution is the proper vertex interaction, which is dominated by the tree-level 
</>i-coupling. The second contribution is the mixing of the fields <p2,s with the 0i. Such a 
mixing of Higgs states occurs at the tree level and is very significant for generic Higgs-boson 
mass matrices, as only a 2% component of 0i will give an effective 02,s-coupling of order 
hf M at t/3 = 50. Since our interest is to assess the significance of the 1-loop singlet-Higgs 
vertex effects, we will mainly focus on variants of the mnSSM and the NMSSM, where the 
mixing of 0i with the other scalars is suppressed. 
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3.1 Decoupling via a Heavy Charged Higgs Boson 

One way to switch off the Higgs-boson self-energy transitions 0i — > 02, s and/or a\ — > a 2 ^s 
is to consider the decoupling of the states 0i and a\ via a heavy charged Higgs boson, with 
M H ± > M S usy- Although the decoupling of 4>i m ay be easily achieved within the MSSM 
for relatively large values of tan/? [16], the situation becomes a bit more involved in its 
singlet extensions. In particular, for the CP-even Higgs sector, one has to arrange that 
both mass-matrix elements (M|) 12 and (M|) 13 vanish, which is more difficult. However, 
this difficulty is not present for the CP-odd Higgs sector, where the mass-matrix element 
(M|) 12 could vanish for certain choices of the model parameters, thereby decoupling the 
CP-odd state a\ from the rest of the mass spectrum. 

3.1.1 Mixing in the mnSSM 

Before we discuss the Higgs-mixing effects, we first give a brief overview of the Higgs sector 
of the mnSSM. The renormalizable part of the mnSSM superpotential is given by 

W mnSS M = hE T d iT 2 LE + h d Hjir 2 QD + h u Q T iT 2 H u U 

+ \SHjir 2 H u + t F S . (3.2) 

In (13. 2p the term linear in S is induced by supergravity quantum effects from Planck- 
suppressed non-renormalizable operators in the Kahler potential and superpotential. De- 
pending on the discrete R symmetries imposed on the theory, the effective tadpole param- 
eter tp and its associate soft SUSY-breaking term ts S may be generated at loop levels 
higher than 5 and can both be of order Msusy- These two interactions are essential to 
break the unwanted PQ symmetry. Further details related to the tadpole generation and 
the breaking of the PQ symmetry may be found in [15]. 

In the mnSSM, the tree-level CP-odd mass matrix Mp reads 

K) 12 = v s ^ CpMl + m?2 ) ' 

(M P ) = t,(s^Ml + m\ 2 ) + ^, (3.3) 

where subscript 1 refers to the CP-odd state a = —s^ai + 0,302 and the subscript 2 to the 
CP-odd state ag- Moreover, M a is the would-be MSSM pseudoscalar Higgs mass, related 
to the charged Higgs-boson mass by M% = M\± — + \ \ 2 v 2 at the tree level, and m\ 2 is 
related to the effective superpotential tadpole tp by m\ 2 = Xtp. It is important to comment 
that the dominant scalar top and bottom corrections to the CP-odd mass matrix Mp can 
all be absorbed to M a and so do not modify its tree-level form. 
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Figure 4: Masses of the Hi (solid line), H 2 (dashed line) and A\ (dot-dashed line) bosons 
in the mnSSM with n = 110 GeV, \t s /n = 150 2 GeV 2 and m\ 2 = -0.5 TeV 2 . The values 
of other soft SUSY -breaking parameters are given in k3. 



By analogy, the tree-level CP-even mass matrix Mj is given by 

(Ml) n = c 2 M 2 + s}Ml, 

(M|) i2 = - spcp {Ml + Mf - AV) , 

( M ^13 = -^(4^ + ^12-2^) > 

(M|) 22 = slMl + clMl, 

(M|) 23 = --^( S/3 c 2 M a 2 + C/3 m? 2 -2 S/3yU 2 ) , 

(M 2 ) 33 = S/3C/3 Q 2 ( S/3C/3 M 2 + m 2 2 ) + ^, (3.4) 

where the subscripts 1, 2, 3 refer to the CP-even states 0i,2,5, respectively. In our numerical 
analysis we also include the 1-loop corrections to M| due to both (s)top and (s)bottom 
loops, which play an important role both for intermediate and large values of tp. 
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Figure 5: The SM-normalized couplings H^bb, H 2 bb and Aibb in the mnSSM, as functions 
of \, for the same model parameters as in Fig. \Q Also shown are the corresponding cou- 
plings Hibb and H 2 bb with the singlet threshold correction A b s = for comparison. The 



pseudoscalar couplings are approximately zero when A^ s 



0. 



To cancel the mixing of state a\ with the other CP-odd Higgs states, we only need to 
choose m\ 2 , such that it satisfies the relation 

spCpM 2 a + m\ 2 = . (3.5) 

This relation also approximately cancels (M|) 13 in the large tp limit. We shall enforce this 
constraint at all times when considering the mnSSM. 

In Fig. H] we display the masses of the two lightest CP-even Higgs bosons Hi and 
H 2 and the lightest CP-odd Higgs boson A\. Since we have taken M H ± = 5 TeV and 
Ats//i = 150 2 GeV 2 , the remaining physical Higgs states, H 3 ~ (pi and A 2 ~ a, are heavy 
of order Mh±- We observe that the lightest Higgs boson mass Mh 1 goes well below the 
LEP limit from direct Higgs searches, for large values of A > 0.3. Hence, the coupling A 
has to be smaller than about 0.3 for a phenomenological viable model. In Fig. we then 
show the dependence of the 6-quark Yukawa couplings 2bb and g^ lbb , for the aforemen- 
tioned scenario. We find that the CP-even Yukawa couplings 2bb receive appreciable 
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contributions from the tree-level mixing of the state with (j) 2i s, which is competitive to 
the loop-induced Yukawa couplings Af 2 ' s . On the other hand, the coupling ~ 9a s bl * s 
completely dominated by the 1-loop contribution A^ s . For moderate values of A ~ 0.3, we 
find that g^ h i ~ 0.15, so the ^i&fr-coupling is ~ 15% of the SM Higgs boson coupling to 
the b quark. Moreover, the decay A\ — > bb is expected to be the dominant decay channel 
in this specific scenario of the mnSSM. 

For completeness, we present numerical estimates for the effective Higgs-boson cou- 
plings to the W ± and Z bosons. These can be determined by the effective Lagrangian 

Chvv = 9u,M w y EgH i w(H i W+W-' lt + ^H.Z.zA , (3.6) 

i=l \ C w / 

where c w = ^Jl — = M w /M z and 

g Hi vv = Ogc? + Olsp. (3.7) 

In the scenario of the mnSSM specified above, the SM-normalized couplings gn x 2 vv are 
shown in Fig. [61 Combining the results of this last figure with Figs. H] and [5l we observe 
that both the lightest CP-even Higgs bosons H\ 2 will predominantly decay into b quarks, 
for small values of A, e.g. for A ~ 0.05. For larger values, i.e. for A > 0.1, only the Hi boson 
will decay into b quarks, whereas the heavier one H 2 will decay into W ± bosons. 



3.1.2 Mixing in the NMSSM 

We now turn our attention to the NMSSM. The superpotential of this model is given by 

Wnmssm = hiHjiT 2 LE + h d Hjir 2 QD + h u Q T ir 2 H u U 

+ \SH T d ir 2 H u + . (3.8) 

The difference between the mnSSM and the NMSSM is that the effective tadpole parameter 
linear in S in the former model is now replaced by an operator cubic in S. In addition to 
the superpotential term |5* 3 , there will be a soft SUSY-breaking operator | A K S 3 , which 
needs be considered as well in the calculation of the Higgs-boson mass matrices. 

In the same weak basis as the one considered for the mnSSM, we first present the CP- 
odd Higgs-boson mass matrix Mp. At the tree level, Mp may be conveniently expressed 
as 

{m 2 p ) u = Ml, 

( M ^) 22 = ^ w (« M « - V 2 ) + ( 3 - 9 ) 
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Figure 6: Effective vector boson couplings for the Higgs bosons Hi (solid line) and 
H 2 (dashed line) in the mnSSM, for the same model parameters as in Fig\4\ 



Here s^c^M 2 = (jaA\ — f/U 2 ) at tree level. In similar fashion, the individual matrix ele- 
ments of the tree-level CP-even Higgs-boson mass matrix Mj are written down as follows: 

(M|) n = cjMf + s^Ml, 

(M|) i2 = - w (m 2 + M 2 -AV) , 



(M|) 22 = s%Ml + c 2 M, 
{Ml) 13 = (slcpMl - 2c^ 2 - ~ S/3/ / 2 ) 

{Ml) 23 = ~ (^c 2 M 2 - 2s m 2 - jc^ 



( M i) 33 = (wM a 2 + ~/) - + 4^/x 2 . (3.10) 

Following the same rationale as in the mnSSM, we suppress the mixing of the CP-odd 
state a ~ — a\ with as by requiring that (Mp) 12 = 0. Hence, we find that 
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Figure 7: The SM-normalized couplings H^bb (solid line), A^bb (dashed line) and 
A\t + t~ (dot-dashed line) in the NMSSM, as functions of A, where the parameters k and 
A K are constrained by A3.11\) and A3.13]) . We also show the coupling H\bb with the singlet 
threshold correction Af" = 0. Both pseudoscalar couplings Aiff are approximately zero 
when Af = 0. 



Note that for positive A, we must have negative k to fulfil the above constraint. Substi- 
tuting 03. Ill) into the expression for (Mp) 22 given in 03. 9I) . we get the mass for the lightest 
CP-odd Higgs boson Ai, i.e. 
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M 2 P 
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(3.12) 



For our benchmark value of \i = 110 GeV and for moderate A ~ 0.6 with tp = 50, the 
positivity condition on M\ gives an upper bound on A K ~ 5 GeV. The simplest option 
would be to set A R = 0. The mass of A\ then strongly depends upon the value of A, through 
the factor v/vs- In order to examine the effect of the threshold corrections on very light 
singlets across a larger range of A, we instead allow A K to vary such that 



v 



spcpix - A K 



0.05 ~ 0.06 



(3.13) 
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This, together with the above constraint (13.111) . gives a singlet pseudoscalar mass of Ma 1 = 
6 ~ 9 GeV for M#± = 2 TeV, across the full range of A. The couplings of Hi to bb and of 
A\ to both bb and t + t~ pairs are shown in Fig. [7J In this variant of the NMSSM, lightest 
CP-even Higgs-boson mass Mh 1 is within the range 120-140 GeV across the full range of A. 

The condition (13.111) can be re- written as a constraint on the parameter A\, 

A X = -2^ = 2 -^Ml (3.14) 

In view of (13.141) it is clear that a singlet dominated light pseudoscalar is present in the 
NMSSM spectrum if Ax ~ 200 GeV and A K ~ 5 GeV. This can be naturally arranged in 
gauge or gaugino mediated SUSY breaking scenarios, where A\ and A K are zero at tree 
level. Quantum contributions from gaugino masses produce non-zero values at the one- and 
two-loop level respectively, leading to the approximate scales shown above if the gaugino 
masses are of the order 100 GeV [11]. Although we have artificially set (Mj,) 12 exactly to 
zero in order to make the effect of the one-loop correction explicit, similar considerations 
apply to the more natural scenario of small but non- vanishing pseudoscalar mixing. 

In the last years, there has been some interest in the phenomenology of light Higgs 
pseudoscalars, which may provide an invisible decay channel for a light SM-like Higgs 
boson. If these CP-odd scalars have a large singlet component, it is possible for them 
to escape experimental bounds [17]. It is clear that the threshold corrections can have a 
significant effect on the branching ratios of a light CP-odd singlet scalar for moderate to 
large values of A. Previous studies have considered detection of these particles through 
decays to photon pairs as the dominant mode [18] in the limit of vanishing singlet-doublet 
pseudoscalar mixing. However, our analysis shows that this need not be the case, and the 
impact of the hadronic decays of A\ in so-called "invisible Higgs" scenarios should still be 
considered even in this limit. 

3.2 Decoupling by Tuning A in the mnSSM 

In the mnSSM the turning off of the singlet-doublet mixing of the CP-odd Higgs scalars 
leads automatically to a suppression of the mass matrix element (M|) 13, which in turn 
implies a small mixing of with Nevertheless, it is also possible to decouple 0i from 
the other CP-even Higgs state 2 by tuning A. Hence one may choose a value for A, such 
that it is (M|) 12 = 0, or equivalently 



where 5 ra( j represents radiative corrections to the tree level mass matrix. Evidently, in the 
absence of any fundamental reason, such a scenario should be considered to be somewhat 
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Figure 8: Effective b-quark Yukawa couplings of the Higgs scalars H\, H 3 and A 2 , as 
functions of \t s /n , where fj, = 110 GeV, M H ± = 120 GeV, m\ 2 = -292 GeV 2 and 
X = 0.7895. Also shown are the corresponding couplings Hibb and H^bb with the singlet 
threshold correction Af s = for comparison. The pseudoscalar couplings are approximately 
zero when A^ s = 0. 



contrived, as it relies upon an unnatural cancellation of different terms to a relatively high 
level of precision. 

In Fig. [8] we present numerical estimates of the Yukawa couplings for this contrived 
model. As input values, we take M#± = 120 GeV and A = 0.7895. The fe-quark Yukawa 
couplings of Hi ~ <p 2 , H 3 ~ 4>s an d A 2 ~ as are plotted as functions of singlet mass 
parameter Xts / /i. The threshold corrections are independent of the singlet mass scale 
\ts/n, so we expect a flat line when these dominate, as is the case for the CP-odd Higgs 
field A 2 . Instead, the respective couplings of the CP-even Higgs bosons show a noticeable 
variation with increasing \ts/n, which originate from the fact that the mass matrix elements 
(M|) 13 and (M|) 23 do not exactly cancel in this contrived model. As a consequence, the 
Yukawa coupling constant gujib ~ 9<t> 2 bb receives an appreciable contribution thanks to 
the mixing between the doublets via an intermediate singlet state, i.e. <p 2 — > <ps — ¥ 0i- 
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Considering LEP constraints, this scenario of the mnSSM is only viable for heavy Higgs 
singlets with masses around 300 GeV. 



4 Conclusions and Future Directions 

Minimal extensions of the MSSM generically include singlet Higgs bosons. Although singlet 
Higgs bosons have no direct or proper couplings to the SM particles, their interaction with 
the observed matter can still be significant as a result of two contributions. The first one is 
their mixing with Higgs doublet states, which is the one often considered in the literature. 
The second contribution is novel and persists even if the Higgs doublet- singlet mixing is 
completely switched off. It results from gluino, chargino and squark quantum effects at the 

1- loop level and has been the focus of this paper. 

In this article we have derived an effective Lagrangian which describes the interac- 
tions between the Higgs bosons and the down-type quarks and leptons in CP-conserving 
singlet extensions of the MSSM. We have found that the loop-induced singlet Higgs-boson 
couplings are enhanced for large values of tan/3 > 40. We have examined the effects of 
these couplings on the phenomenology of two such models, the mnSSM and the NMSSM. 
Although the mixing of <f>\ with the other CP-even Higgs scalars <p2,s is the leading effect for 
most of the parameter space, we have found that the threshold corrections to the Yukawa 
couplings remain relevant and can play an important role in the phenomenology of possibly 
light singlet Higgs scalars. 

In the absence of a Higgs doublet-singlet mixing, the 1-loop quantum effects we have 
been studying here will be the only means by which the CP-odd singlet may couple to 
quarks and leptons. For a sufficiently light CP-odd singlet scalar, with a mass below 
the squark threshold, the loop-induced Yukawa couplings will provide its dominant decay 
channel into b quarks. This has important phenomenological implications, since a SM-like 
Higgs boson will no longer decay invisibly into a pair of light singlet pseudoscalars A\A\. 
The threshold corrections that we calculated here will give rise to the leading decay mode 
A\ — > bb. In fact, this decay channel depends on the mass of Ax, and can have a coupling 
strength of up to ~ 1/2 of the corresponding SM coupling, i.e. g^ lbb ~ 0.5. 

There are several possible new directions for future study. For instance, one may 
lift the assumption of CP conservation in the singlet extensions of the MSSM. Then, light 
CP-odd scalars may contribute to electron and neutron electric dipole moments at the 

2- loop level [19]. It would be interesting to study the impact of those contributions in 
the presence of CP-violating threshold corrections. Another possible direction for future 
investigations will be to calculate the off-diagonal couplings of the singlet Higgs bosons 
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to down- type quarks [20]. Our effective Lagrangian presented here may be generalized to 
include these flavour- changing neutral-current (FCNC) interactions of the Higgs bosons to 
quarks. It would be particularly valuable to explore the impact of the singlet Higgs-boson 
FCNC effects on K- and 5-meson observables. 
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A Background-Field-Dependent Masses and Mixing 
Angles 



In this appendix we present our notation and conventions for the masses and mixing angles 
which enter the calculations of Sections 12.11 and 12.21 We have divided the appendix into 
two subsections. The first gives the relevant squark mass parameters, whilst the second 
contains the corresponding information for the chargino sector. 

A.l Squark Masses and Mixing Angles 

Neglecting D-term contributions, the scalar top and bottom mass matrices may be written 
in the (qx, gi?) T basis as 



Ml 



Mq + h 2 t \$° 2 \ 2 h t A t $f - h t \S§°* 



h t A t §° 2 - htXS*^ M 2 + h 2 \^? 



012 
2 1 



(a.i; 



Ml 



(A.2) 



M§ + /ig|$o|2 h b A h <&l - h b \S<S>° 2 
h h A b $f -h h \S*<5>° 2 * Ml + h\\§>\\ 2 

Where the fields <3>° and 5* have been defined after (12. ip . These matrices are diagonalized 
by unitary matrices which may be parameterised as 



where 




and c„ = cos 9„ 



cos 9 q sin 6 q 
— sin 9 q cos 6 q 



s t = sin 9 t 



s b = sin 9 b 




S q = arg((^)J 



\ 



m| - M 2 - hi m\ 2 



m 



t-2 



ill} 



"i - M Q - fe wi 2 

1 2 

m{ — mf 

02 »1 



The field-dependent scalar quark masses are given by 



(A.3) 



(A.4) 



5,2 



m 



Cl(2) 



^2 



m 



'1(2) 



M 2 + M 2 + 2h 2 \^° 2 \ 2 + 



M 2 - M 2 ) + 4/if \A t $° 2 - AS*$?| 2 



M 2 + M 2 h + 2h 2 M 2 + {-)\J(M 2 Q - M 2 ) + Ah 2 b \A b ^ - A5$ 



012 
2 1 



CA.5) 
.(A.6) 
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A. 2 Chargino Masses and Mixing Angles 



The Lagrangian describing the chargino masses and their Yukawa interactions to Higgs 
bosons is given by 



-iw L ,h 1L ) M 



-iw£ 
ht L 



h.c, 



M 



with 

M 2 g w <$>° 2 * 

The mass matrix M may be diagonalized by the bi-unitary transformation 



(A.7) 



(A.9) 



where Md = diag(m Xl , m X2 ). Both the chargino masses and the elements of the unitary 
matrices IA and V depend explicitly on the Higgs background fields $5 2 an d S. 

We may parameterise the unitary matrices U, V as 



V 



cos6*~ 
sin 9~ 



— sin 6 

cos 9' 



cos 9 + — sin # H 
sin# + cos9 + 



1 

e iS ~ 



e l 



(A.10) 
(A.H) 



with 



6* 



arg (M 2 $?* + AS*$°* 
arg (M 2 $1* + AS*?* 



sin 



\ 



sin 



0+ 



m xi 


-Ml 




m\ 2 




m L 


-m 2 Xl 




m 2 Xl 


-M| 


- V 


$0|2 






-m 2 Xl 





(A.12) 



and = cos^ = J 1 — (s^ 2 . The chargino masses are given by 
± [M| + A 2 |5| 2 + 2gl (|$?| 2 + |$^| 



X 2 (l) 



(A.13) 



(M| + A 2 |S| 2 + 2g 2 w (|$?| 2 + |$Q| 2 )) 2 - 4|M 2 AS - 2^$?$°*| 2 
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